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ABSTRACT 

We present new observations of the Kuiper belt analogues around HD 38858 and 
HD 20794, hosts of super-Earth mass planets within 1 au. As two of the four nearby 
G-type stars (with HD 69830 and 61 Vir) that form the basis of a possible correlation 
between low-mass planets and debris disc brightness, these systems are of particular 
interest. The disc around HD 38858 is well resolved with Herschel and we constrain 
the disc geometry and radial structure. We also present a probable JCMT sub-mm con¬ 
tinuum detection of the disc and a CO J=2-l upper limit. The disc around HD 20794 
is much fainter and appears marginally resolved with Herschel , and is constrained to 
be less extended than the discs around 61 Vir and HD 38858. We also set limits on the 
radial location of hot dust recently detected around HD 20794 with near-IR interferom¬ 
etry. We present HARPS upper limits on unseen planets in these four systems, ruling out 
additional super-Earths within a few au, and Saturn-mass planets within 10 au. We con¬ 
sider the disc structure in the three systems with Kuiper belt analogues (HD 69830 has 
only a warm dust detection), concluding that 61 Vir and HD 38858 have greater radial 
disc extent than HD 20794. We speculate that the greater width is related to the greater 
minimum planet masses (10-20 M 0 vs. 3-5 M®), arising from an eccentric planetesi- 
mal population analogous to the Solar System’s scattered disc. We discuss alternative 
scenarios and possible means to distinguish among them. 

Key words: planetary systems: formation — planet-disc interactions — circumstellar 
matter — stars: individual: HD 20794 — stars: individual: 61 Vir — stars: individual: 
HD 69830 — stars: individual: HD 38858 


1 INTRODUCTION 

The first evidence for the existence of a planetary system other 
than our own was arguably the image of the debris disc around 
/3 Pictoris JSmith 8c Terrilel IT984h . With the contemporaneous 
discovery of many other examples of th(V‘Vega phenomenon” 
with the IRAS satellite dAumann et al.lll984h . it was soon re¬ 
alised that the dust in these system s was the product of colli - 
sions between larger planetesimals feackman & Parescell 19931) . 
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the same building blocks thought to build the planets now seen 
around hundreds of stars. 


Despite a clear conceptual connection between planets and 
debris discs, both being a product of some degree of planet 
formation, there has until recently been surprisingly little ev- 
idenc e of statistically significant correlations between the two 
(e.g. Brvden et al.ll2009l: iKospal et al.l[2009l : iKennedv & Wvattl 
12012 : iMorales et al J 1201 2h ~ Ho wever, with many ne arby stars 
observed with both the Spitzer dWerner et al.l 120041) and Her- 
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,k:V?<?{ 3 dPilbratt et al .1120 1 (j) observatories, and the continual in- 
crease in the sensitivity of radial velocity su rveys for p lanets, 
the first correlations have begun to eme rge dWvatt et al.ll2012l : 
iMarshall et al.ll2014h . lWvatt et al.1 d2012h reported that of 6 out 
of the 60 nearest G-type stars with close-in low-mass planets 
and no known planets above a Saturn-mass, 4 have relatively 
bright debris discs (i.e. at least an order of magnitude brighter 
than our Edgeworth-Kuiper belt; HD 20794, 61 Vir, HD 69830, 
HD 38858). This high detection rate is in contrast with a de¬ 
tection rate of ~15% around normal Sun-like main-sequence 
stars dTrilling et al.1 120081) . which is comparable with that for 
stars with more massive planets, suggesting there is some kind 
of link between low-mass planets and the brightness of the de¬ 
bris in those systems. A decreased ability of low-mass planets 
to eject planete simals compared to gas giants is a possible ex¬ 
planation (e.g. lPavne et ~ail2009i: I Wyatt et al.ll2012h . 

With only a handful of low-mass planet + disc systems, the 
cause of such a correlation remains unclear, so characterisation 
of the four G-type systems that contribute to this correlation is 
a first step towards an understanding. Of those four, resolved 
imaging of 61 Vir was presented and discussed bv IWvatt et al.l 
J2012h . and HD 69830, which only_hosts warm dust that was not 
detected by Herschel (Marshall et al. 201 4), has been discussed 
in numerous stud ies (Beichman et al. 20 051 : lEisse et al .1 [20071 : 
ISmith et aPl2009l : iBeichman et al.|[201~lh . Here, we present new 
Herschel PACS observations of the remaining two systems; 
HD 38858 and HD 20794. The former hosts a bright and well 
resolved debris disc, so was also observed with SCUBA-2 and 
RxA3 on the JCMT. The latter hosts one of the faintest debris 
discs known, and at only 6 pc and with three low-mass planets 
is clearly a system of high interest. In what follows we present 
and describe the observations, derive constraints on the disc spa¬ 
tial structure, present upper limits on undiscovered planets us¬ 
ing radial velocity (RV) data, and discuss the disc structure and 
possible relations with the planets. 


2 THE HD 38858 AND HD 20794 SYSTEMS 

We begin by describing the previously known stellar and plan¬ 
etary properties of our systems. HD 38858 (HIP 27435) is a 
nearby (15.2 pc) Sun-like (G4V)_star. Stellar age estimates vary 
from as young as 200 Myr (Casagran de et al.11201 ll), to as old as 
2.32 to 8.08 Gyr (iTakeda et al .11 20071) . and 9.3 Gyr dSousa et al.1 
2010). The star is known t o host a cool debris di sc, which was 
first detected with Spitzer feeichman et al.ll2006i) . As noted by 
iKrist et ail (120121) . the disc was found to be marginally resolved 
at 70 pm with Spitzer with an inclination of 48°, position angle 
of 56°, and a radius 135 au (9", no uncertainties given). The star 
was observed with the Precision Int egrated Optics Near Infrared 
ExpeRiment (PIONIER. lLe Bouauin et ai1l201 ll) in search of a 
resolved excess, and the 3 a upper limit on the disc/star flux ra¬ 
tio in the H-band is 1% fertel et alJl2014h The j>tar has been 
reported to host a super-Earth mass planet dMavor et al.11201 ilk 
for which we provide an update on the orbit below (Marmier et 
al. in preparation). 

HD 20794 (HIP 15510) is a nearby (6 pc) Sun-like (G8V) 


1 Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 


Table 1. Herschel PACS observations of HD 38858 and HD 20794. 
All PACS observations use either the 70 or 100 pm bands, and always 
include the 160 pm band. 


Target 

Obslds 

Date 

A (pm) 

Duration (s) 

HD 20794 

1342216456/7 

20 Mar 2011 

100 

2 x 1686 

HD 20794 

1342234096/7 

14 Dec 2011 

70 

2 x 445 

HD 38858 

1342242537/8 

28 Mar 2012 

70 

2 x 1686 


star. Stellar age estim ates vary from as young as 3.9 Gyr 
(ICasagrande etal.1201 ll) . through inter mediate ages of 5.2 to 6. 2 
Gvr d Vicanl 120 1 2h and up to 11.3 Gyr (iHolmberg et al.ll2009l) . 
The star has been found to host two robust low-mass planets, 
with a possible third planet JPepe et al.l201lh . The star was sub¬ 
sequently found to host a debris disc using archival Spitzer data, 
as described bv lWvatt et al. ( 2012 ). who also reported a detec¬ 
tion with Herschel that we elaborate on below. This star was 
also observed with PIONIER, resultin g in a s ignificant H -band 
disc/star flux ratio of 1.64 =t 0.37 fertel et al .11 20141) . 


3 DEBRIS DISC OBSERVATIONS 

Here we present new observations of HD 38858 with sev¬ 
eral different instruments. It was observed by Herschel 
feilbratt et al.l l2010l) under the auspices of the Search for 
Kuiper-belts ARound Planet-host Stars (SKARPS) survey (e.g. 
iKennedv et al.ll2013l) . and by the new Submillimetre Common 
User Bolometer Array (SCUBA-2) as part of the SCUBA-2 
Observations of Nearby Stars (SONS) James Clerk Maxwell 
Telescope (JCMT) Legacy survey feanic et al .1120131) . Finally, 
HD 38858 has also been observed in search of CO J=2-l emis¬ 
sion using the RxA3 receiver on the JCMT. 

We also present a detailed analysis of observations of 
HD 20794 with the Herschel PACS instrument, observed as 
part of the Disc Emission via a Bias-free Reconnaissance in the 
Infrared/Submillimetre (DEBRIS) survey. The disc around this 
star is very faint compared to HD 38858, so it was not selected 
for the any of our JCMT programmes due to the high probability 
of a non-detection. 


3.1 HD 38858 

3.1.1 Herschel PACS 

HD 38858 was observed with the H ersch el Photodetector Array 
Camera and Spectrometer (PACS. lPoglitsch et al.ll2010l) instru¬ 
ment at 70 and 160 pm as part of the SKARPS open time pro¬ 
gramme (see TableQ}. The 70 pm image of HD 38858 is shown 
in Fig. Q] The left panel shows the original image and the stel¬ 
lar emission has been subtracted from the right panel. Most of 
the 70 pm emission clearly originates in the disc, and the full- 
width at half-maximum (FWHM) PACS beam size at 70 pm 
is about 6" so the disc is clearly resolved. The disc is approx¬ 
imately 15" across, so the system distance of 15.2 pc implies 
a disc diameter of around 225 au, similar to the size estimated 
from Spitzer IKrist et al.ll2012l) . We return to the disc properties 
and geometry in section I47TI 

Measuring the brightness of the star+disc in the PACS im¬ 
age is fairly straightforward for HD 38858. The disc is well re¬ 
solved so we use a 20" radius aperture that is centred on the 
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Figure 1 . Herschel PACS 70 /xm images of HD 38858. The left panel shows the original image and the star has been subtracted from the right panel. 
The star location is indicated by the white cross. Contours are at 3, 5, 10, and 25 times the pixel RMS level of 2 x 10 -2 mJy arcsec -2 (70 gm pixels 
are 1 arcsec 2 ). The PACS beam size at 70 gm is about 6 ", so the disc is well resolved. 



peak emission location (and about 3" from the expected stel¬ 
lar position, so consistent with the pointing uncertainty). After 
aperture correction we find a flux density of 190 zb 6 mJy at 70 
gm, and of 141 ± 20 mJy at 160 gm. The uncertainty is esti¬ 
mated by the mean value of a number of 20" apertures placed 
randomly in high coverage regions in each map, and is added 
in quadrature to the calibration uncertainty (which is included 
above, and dominates the 70 gm uncertainty). The uncertainty 
is several times higher at 160 gm than expected simply from the 
PACS integration time due to a high Galactic background level. 
This high background level is expected in the Orion Complex, 
and HD 38858 was excluded from the DEBRIS survey based on 
such an expectation ([Phillips et al.ll20ld) . These measurements 
are in excellent agreement withji 70 gm Spitzer measurement 
of 190 =b 9 mJy dLawler et al .1120091) . 


3.1.2 JCMT SCUBA-2 continuum imaging 

As a star with a disc bright enough for a possible detection at 
sub-mm wavelengths, HD 38858 was i ncluded in the SONS sur¬ 
vey, which uses SCUBA-2 ( Holland et al. 2 013b to detect and 
image discs at 450 and 850 gm. HD 38858 was observed for a 
total of 10 hours, made up of 20 half-hour integrations. We used 
the constant speed DAISY pattern, which provides uniform ex¬ 
posure time coverage i n the central 3 arcmin-diameter region of 
a field llHolland et al.1120131) . 

The SCUBA-2 data were reduced using the Dynamic 
Iterative Map-Maker within the STARLINK SMURF pack¬ 
age dCha pin et al. 20131) c alled from the ORAC-DR automated 
pipeline ( Cavanagh et al.lE008l) . The map maker used a config¬ 
uration file optimized for known position, compact sources. It 
adopts the technique of “zero-masking” in which the map is 
constrained to a mean value of zero (in this case outside a ra¬ 
dius of 60" from the cente r of the fi eld), for all but the final 
interation of the map maker dChapin et al.120131) . The technique 
not only helps convergence in the iterative part of the map¬ 
making process but suppresses the large-scale ripples that can 


produce ringing artefacts. The data are also high-pass filtered at 
1 Hz, corresponding to a spatial cut-off of ^150" for a typical 
DAISY scanning speed of 155"/s. The filtering removes resid¬ 
ual low-frequency (large spatial scale) noise and, along with 
the “zero-masking” technique, produces flat and uniform final 
i mage s largely devoid of gradients and artefacts dChapin et al.1 

To account for the attenuation of the signal as a result of 
the time series filtering, the pipeline re-makes each map with 
a fake 10 Jy Gaussian added to the raw data, but offset from 
the nominal map centre by 30" to avoid contamination with any 
detected source. The amplitude of the Gaussian in the output 
map gives the signal attenuation, and this correction is applied 
along with the flux conversion factor derived from the calibrator 
observations. The final images were made by coadding the 20 
maps using inverse-variance weighting, re-gridded with 1" pix¬ 
els at both wavelengths. The final images at both wavelengths 
have been smoothed with a 7" FWHM Gaussian to improve 
the signal-to-noise ratio. The FWHMs of the primary beam are 
7.9" and 13.0" at 450 and 850 /im, respectively. 

Fig. m shows the SCUBA-2 images overlaid as contours 
on the PACS 160 gm data. At 450 /xm a 4 a significant peak 
is seen about 10" to the North of HD 38858. This detection 
is consistent with being unresolved, and has a flux density of 
49 zb 11 mJy. The relatively low S/N means that the detected flux 
may actually originate closer to HD 38858 than 10", but given 
the known high background level it is likely that this detection is 
not associated with HD 38858 and we do not consider it further. 

At 850 gm the SCUBA-2 data show a clear detection. The 
peak location (6a) is offset from the star location by about 
10" to the East, and the centre of the 5 a contour offset by about 
5". A background source is seen to the South, and probably 
corresponds to emission detected at 160 /xm at about the same 
location with PACS. The co-location of these sources suggests 
that the accuracy of the JCMT pointing for these observations 
is within a few seconds of arc, and splitting the SCUBA-2 data 
in half yields two images with similar structure. More broadly, 
we find no systematic pointing errors in SONS observations, for 
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Figure 2. SCUBA-2 images of HD 38858 at 450 (left panel) and 850 gm 
contours (solid are positive, dashed are negative), and the rainbow colourec 
location is indicated by the white cross. As indicated by the circles in the 
and 850 fi m. 
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(right panel). The contours show the SCUBA-2 images with 3, 4, 5 and 6cr 
1 background image in both panels shows the 160 gm PACS image. The star 
lower right of each panel, the SCUBA-2 beam FWHM is 8 and 13" at 450 


example the pre- and post-observation pointing calibrations are 
within about 3" across the entire SONS survey. The 850 fim flux 
density at the peak position is 4.9 =t 0.7 mJy/beam, and at the 
stellar position 3.7 ±0.7 mJy/beam. The integrated flux within 
30 // radius of the peak location is 11.5 ± 1.3 mJy, the larger flux 
indicating that the emission covers a few JCMT resolution ele¬ 
ments. The 3-5cr contours cover the location of HD 38858, so 
the image could include both Orion Complex background emis¬ 
sion to the East of the star and disc emission centred on the star. 
If the contamination is not severe, the disc flux lies between 3.7 
and 11.5 mJy. 

An alternative possibility is extra-galactic confusion; 
given about 10 3 galaxies degree -2 brighter than 5 mJy (e.g. 
iBlain et alii 19991) . there is a 2% chance that a source as bright 
as the peak would appear within an aperture 15" in radius cen¬ 
tred on HD 38858, and a much smaller chance that such a source 
would be well centred on the star. Thus, the star-centred emis¬ 
sion could arise from a galaxy with ^0.2% probability, or the 
offset peak could arise from a galaxy with 2% probability. In 
the latter case the star-centred emission still requires an expla¬ 
nation, and adding a second galaxy lowers the probability sig¬ 
nificantly. Thus, in either case it is unlikely that the emission lo¬ 
cated at the stellar position originates from a background galaxy. 

Given the small chance of contamination from background 
galaxies, and that some contribution from the Orion Complex is 
likely, we consider the SCUBA-2 emission a probable detection 
of emission associated with HD 38858. Our favoured interpreta¬ 
tion is that because the 3-5cr contours cover the stellar location, 
850/xm emission from the disc is indeed detected with a flux 
of at least 3.5 mJy (i.e. >5x the point source uncertainty) and 
that the emission is shifted to the East by extended background 
contamination that also includes the source to the South. We il¬ 
lustrate this interpretation further after constructing models of 
the HD 38858 disc below in section l4X3l where we find an 850 
m disc flux of 7.5 ± 2 mJy. 



Figure 3. JCMT RxA3 spectrum for HD 38858. The velocity of 
HD 38858 is 31.54 ± 0.04 km s _1 . The dark red line shows the veloc¬ 
ity profile for a ring of material at 200 au at the detection limit of 1 Jy 
km s -1 . 


3.1.3 JCMT CO J-2-1 line emission 

The RxA3 receiver on the JCMT was used to search for CO 
J=2-l (230.538 GHz) emission towards HD 38858. The grid- 
chop mode was used to stare at the target while beam-switching 
on the sky. Standard calibration was done using hot, cold and 
sky loads, and inspection of spectra of standard sources. The 
co-added spectrum consists of 27 observations of 20 minutes 
duration each, over 5 nights between 12 Nov 2013 and 24 Dec 
2013. The correction ? 7 (main-beam) (T m b = /^mb) at the 

time was between 0.5 and 0.6 from measurements on Mars. The 
22 individual outputs of the ORACDR pipeline were co-added 
using the ’wcsmosaic’ task in KAPPA. The output data have 
units of TX in K versus heliocentric velocity (radio definition) 
in km s -1 , and are converted to Jy with multiplication by a fac- 
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tor of 27.4 Jy K -1 . The spectral ’cube’ (1x1x2048 points) was 
exported to Starlink SPLAT for further reduction. The spectrum 
was rebinned by a factor of 5 to a 0.2 km s -1 channel width, and 
a low-order polynomial fit (that excluded the absorption feature 
described below) subtracted. 

The RxA3 spectrum is shown in Fig. [3] in Jy channel -1 . 
A spectral line associated with the star should be centred 
aroun d the heliocentric stellar radial velocity of 31.54 km s -1 
dNidever et all 2002b . To guide the eye the dark red line shows 
an example line profile with an integrated flux of 1 Jy km s -1 , 
calculated as the histogram of velocities for material uniformly 
distributed around a circular narrow ring at 200 au, with the in¬ 
clination derived from the 70 gm PACS image. There is no clear 
line emission at the expected velocity of HD 38858, but a large 
negative feature at 26 km s -1 , which likely originates from the 
Orion Complex background, where CO J= 1-0 emission is seen 
at 23.5 to 30 km s -1 (heliocentric. !Wilson et all 20051) . This dip 
in the spectrum makes estimation of the baseline difficult; the 
flux in the range 27-30 km s -1 may suffer from absorption and 
thus the flux near the velocity of HD 38858 could be higher than 
estimated in Fig. [3] Future observations that find the origin of 
this line and subtract it will help resolve this issue. 

We estimate the uncertainty on the measurement by fitting 
line profiles like that shown in Fig. [3] for a range of velocity 
widths. We first randomized the bins in the spectrum multiple 
times, fitting lines of a range of widths each time, resulting in 
a standard deviation of about 1 Jy km s -1 for radii from 50- 
300 au. We repeated this exercise, but at a range of locations in 
the non-randomized spectrum, finding similar results. We revisit 
the implications of this upper limit in terms of the gas mass in 
section [4X4] 


3.2 HD 20794 

We now turn to the observations of HD 20794. The disc around 
this star is much fainter than that around HD 38858 so was only 
observed with Herschel PACS (see TableQ}. HD 20794 was first 
observed at 100 gm and 160/im, and following the possible res¬ 
olution of the disc at low S/N, was re-observed at 70 and 160 
/im. 

In contrast to HD 38858, most of the emission observed 
near HD 20794 with PACS comes from the star itself. The 70, 
100, and 160 gm images of HD 20794 are shown in Fig.|4] The 
top row of panels show the 70 gm image, where the star has 
been subtracted in the right panel. The middle row of panels is 
the same, but shows the 100 gm image. The bottom row shows 
the two 160 /im images, the image in the left panel was taken at 
the same time at the 100 /im image, the right panel at the same 
time as the 70 gm image. The star has not been subtracted from 
these bottom panels. 

It is clear from the 70 and 100 gm star-subtracted panels 
in Fig.[4](i.e. top right and middle right) that removing the stel¬ 
lar flux density does not account for all of the emission located 
at the stellar position. There is also evidence that the residual 
emission at the star location is resolved, as the region enclosed 
by the 3 a contour in the star-subtracted 100 gm image is ellip¬ 
tical with a position angle pointing slightly East of North. We 
return to the issue of whether the emission is resolved when 
modelling the disc below in section l4~2l 

The 100 and 160 gm images also show that there is signif¬ 
icant emission that lies to the NE of the star, as marked by the 


“x” symbol. HD 20794 has a high proper motion of 3.0" yr -1 
Eastwards and 0.7" yr -1 Northwards (Ivan Leeuwenl120070 . so 
the nearly year long interval between the 100 and 70 /im obser¬ 
vations means that a background object would move Westward 
(right) relative to the star by about 2.3" between epochs. This 
movement is not easily discernible for the NE background ob¬ 
ject in these images due to a low S/N, but is possibly evident 
from both the change in separation and position angle of the 
stellar and background source emission in the 160 /im images. 

To measure the flux density of the star+disc emission from 
HD 20794 we use point-spread function (PSF) fitting. Though 
the disc may be marginally resolved, it is also confused to some 
degree by the source 10" away to the NE. This method gives a 
measurement that is more robust to confusion, but that will un¬ 
derestimate the flux density slightly if the disc is resolved. For 
the PSF fitting we find 100 =b 5 mJy at 70 gm, 55 =b 3 mJy at 
100 gm, and 27 zb 3 mJy at 160 gm, including calibration un¬ 
certainties. Previous measu rements at 70 /im using MIPS found 
94 zb 14 and 107 zb 4 mJy (iBeichman et al. 1120061 ; iGasnar et alJ 
horn consistent with our measur ement. Our values are also 
consistent with those presented in Wvat t et al.1 d2012h . We re¬ 
visit the derived fluxes when modelling the disc in section FOl 


4 DISC MODELS 

We now consider the properties of the debris discs around 
HD 38858 and HD 20794. We have different levels of infor¬ 
mation about each, with the HD 38858 disc well resolved and 
bright, but the HD 20794 disc faint and only marginally re¬ 
solved at best, so the modelling has a level of complexity that 
reflects this. We first consider the implications of the flux distri¬ 
butions inferred from the star+disc photometry, and then move 
onto modelling of the resolved images. For HD 38858 we also 
use the disc models to interpret the upper limit on CO emission 
and the mass of gas present. 


4.1 HD 38858 

4.1.1 SED 


We begin the modelling of HD 38858 with the flux den¬ 
sity distribution, which we abbreviate to “SED” for simplic¬ 
ity. The SED from optical to sub-mm wavelengths is shown 
in Fig. \5\ The short-wavelength photometry comes from vari- 
ous sources, and includes “heritage” Stro mgren and UBV bands 
(lHauck & Mermilliodl[ l998: Mermilli o J2006l) and space-based 
Tycho bands jHpg et al.1 2000h . but the quality of the pho- 
tospheric prediction ultimately relies onjnid-IR photome¬ 
try from IRAS, AKARI, and WISE (Neugebauer et al] 1 1984 
llshihara et al.1 l20ld : IWright et~akl l2010l) . HD 38858 was ob¬ 
served with the IRS sp ectrograph, and we use t he spectrum from 
the CASSIS database iLebouteiller et al.l 201 ll) . We also include 
Spitzer MIPS photometry from the aforementioned sources. 

To construct an SED model, we fitted photometry short- 
ward of 15 /im with an AM ES-Cond stellar photosphere model 
(iBrott & Hauschi ldt 2005), finding a stellar effective tempera¬ 
ture of 5780 =b 20K and luminosity of 0.83 L©. With only the 
5-35 /xm IRS spectrum, and far-IR photometry at 70 and 160 
/im, the spectral coverage for the disc is poor, which limits our 
ability to constrain the disc spectrum. For the disc we therefore 
simply subtract the expected stellar flux in each band, and fit the 
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Figure 4. Herschel PACS images of HD 20794, North is up and East is left. Contours are at 3, 5, 10, and 25 times the pixel RMS level for each 
observation, though the higher levels are not visible in all images. The surface brightness in the scale bars applies to both images in each row. Top row: 
at 70 gm, the left panel shows the image and the star has been subtracted from the right panel. Middle row: at 100 gm, the left panel shows the image 
and the star has been subtracted from the right panel. Bottom row: at 160 gm, the left panel shows the image taken with the 100 pm image, and the 
right panel shows the image taken with the 70 pm image (neither is star-subtracted). The “+” symbols shows the star position in each panel, and the 
“x” symbols shows the location of a background object. The high Eastward proper motion means that the background object is shifted about 2" to the 
West (right) in the later 70/160 pm images relative to the earlier 100/160 pm images. 


remaining fluxes with a blackbody disc model@ The tempera¬ 
ture of this disc model is well constrained at 60 d= 2K. 

2 The AMES-Cond model spectra are computed to ~lmm, so the 
few percent decrease in the stellar 100 pm flux density relative to a 


The long-wavelength part of the disc spectrum depends on 
our interpretation of the SCUBA-2 observations (see section 

Rayleigh-Jeans extrapolation from the mid-IR is included in our photo- 
spheric prediction. 
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Figure 5. SED for HD 38858. Dots are fluxes and triangles 3cr upper 
limits. Black symbols are measured fluxes and brown symbols are star- 
subtracted fluxes (i.e. disc fluxes, which may cover the measured pho¬ 
tometry). Grey dots are disc fluxes that are consistent with zero, and 
grey triangles 3 a upper limits on these fluxes. The 5780K stellar pho¬ 
tosphere model is shown in blue, the 60K blackbody disc model in red. 
The disc model is a modified blackbody, which assumes a SCUBA-2 
850 gm disc flux of 7.5 ± 2 mJy. 


13.1.21) . Assuming that the 450 /tm emission is not associated 
with HD 38858, and that the 850 gm emission is partly contam¬ 
inated, we find a disc flux of 7.5 mJy, and the disc spectrum is 
near to a pure blackbody, as shown in Fig. \5\ We use a “mod¬ 
ified” blackbody, where we multiply the spectrum by (Ao/A)^ 
beyond Ao gm, with Ao = 210 gm and f3 = 0.25. With only a 
single measurement beyond 160 gm, Ao and /3 are degenerate, 
and their accuracy is also limited by the systematic uncertainty 
introduced by our interpretation of the SCUBA-2 850 gm im¬ 
age. 

Useful disc properties that can be derived from the SED 
are the disc radius assuming blackbody grains, which yields 21 
au (about 3"), and the disc fractional luminosity / = Tdisc/T*, 
which is 8.8 x 10 -5 . As can already be seen from Fig.Q] the 
disc is much larger than this estimate, indicating that the grains 
in the disc are hotter than blackbodies for a given location, as 
expected for grains that emit poorly at long wavelengths. Thus, 
the requirement of a nearly pure blackbody spectrum to account 
for the SCUBA-2 flux is surprising, since if the disc material lay 
at a single radius a blackbody spectrum implies that the black¬ 
body and resolved disc sizes would agree. 

However, blackbody disc spectra have been seen for other 
discs that appear at larger radii that their blackbody tempera¬ 
ture suggest. In the case of AU Mic, the likely explanation is 
that the sub-mm observations detect the parent body distribu¬ 
tion, which may have blackbody properties, but that far-IR ob¬ 
servations detect a halo of smaller grains that reside on high 
eccent ricity orbits due to radiation and/or stellar w ind forces 
(e.g. IStrubbe & Chian3l2006l : iFitzgerald et alJl2007l Matthews 
et al. submitted). Therefore, a discrepancy between the black¬ 
body and observed disc sizes does not warrant the exclusion of 
sub-mm fluxes that result in near blackbody disc spectra. 

Using the disc temperature and the SCUBA-2 observa¬ 
tions, we can also derive the dust mass in the disc. Assuming 
an 850 gm opacity of 45 au 2 M ^ 1 and flux of 7.5 mJy the 
dust mass is 0.008 M®. Because most bodies are expected to 


be much larger than the detected dust, the total disc mass is 
much larger, depending on the size D c (in km) of the largest 
planetesimals. Using equation (15) of Wvattl (2008) and assum¬ 
ing a disc radius of 60 au and a minimum grain size of 1 gm 
yields a total disc mass M to t ~ D°' 5 M®, so about an Earth 
mass for 1 km sized planetesimals. By assuming a stellar age 
and that the debris disc has been evolving for all of this time, 
D c can be estimated, for example using equation (16) of lWvattl 
( 2008b , where the important proportionality here is D c oc t 2 ge . 
If we assume an age for HD 38858 of 1 Gyr and the planetes- 
imal properties found for Sun-like stars bv [Rains et al] ( 2011 ). 
the maximum size of objects that must be present to replen¬ 
ish the disc over the star’s lifetime is 0.04 km, though this size 
is very uncertain because it depends on the square of the stel¬ 
lar age and planetesimal strength. For weaker planetesimals as 
found for A-type stars by IWvatt et all (l2007bl) the size is of or¬ 
der 10 km. Thus, the corresponding disc mass is between about 
0.1 to 3 Earth masses, and only depends linearly on the age and 
object strength. We can also estimate the rate at which mass is 
being lost from the disc as the total mass divided by the age, 
which yields 10 15 to 10 16 kg yr -1 . Because the derived mass 
depends on the age (via D c ), the mass loss rate is independent 
of the assumed age and the main contributor to the uncertainty 
as calculated here is the planetesimal strength. This mass loss 
rate will be used when considering the mass of CO gas below in 
section [5X4] 


4.1.2 Resolved models 

We now construct resolved models of the PACS images of 
HD 38858 to constrain the disc spatial structure. For a given 
disc structure, we first generate high resolution images at each 
observed wavelength, and then convolve these with the tele¬ 
scope response to a point source (i.e. the PSF, which is a cal¬ 
ibration observation of 7 Dra) for comparison with the observa¬ 
tions. All of our Herschel images are created with North point¬ 
ing upwards, so the PSF observations, which have a different 
telescope orientation, must be rotated to the correct orientation. 

The disc is well resolved and our unsuccessful attempts to 
model the 70 gm image with a single narrow ring find that the 
disc has significant radial extension. We therefore use two sim¬ 
ilar models, one in which the disc is simply made up of two 
narrow dust belts, and one where the disc is radially extended. 
Common parameters for both models are the disc inclination 
i and position angle D (East of North), and the disc tempera¬ 
ture profile as defined by Tauc = /tTbb- The factor /t ac¬ 
counts for how much hotter the dust appears compared to a 
blackbody at the same distance from the star, where Tbb = 
278.3 /y/r. That is, fr reconciles the disc spatial structure 
with the SED. All models are consistent with /t = 1.9, similar 
t o the values found for the disc around the Sun-like star 61 Vir 
(IWvatt et ahll2012l) . We use modified blackbody parameters for 
the disc spectrum as described above, though here our goal is to 
model the disc structure as seen by PACS so they are relatively 
unimportant. 

For the two ring model we use belts that are 10 au wide 
with constant optical depths, and the best fitting model is shown 
in Fig. [ 6 ] The disc inclination and position angle are 44° and 
67°, with estimated uncertainties of 5°. These measurements 
are consistent with the geometry reported bv lKrist etai . ( 2012 ) 
from lower resolution Spitzer observations. The two rings are 
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Figure 6. Two ring model of the HD 38858 disc. The columns show the PACS image, convolved model, high resolution model, and residuals from left 
to right. The top row is at 70 gm, and the bottom row 160 jjm. The model has two discrete rings that are both 10 au wide, and are centred at 55 and 
130 au. This model fits the PACS data very well, with the only significant residuals being due to the high background level at 160 gm. 



Figure 7. Two alternative models of the HD 38858 disc. Each square of four panels shows the high resolution model and residuals at 70 and 160 gm. 
The left set of panels shows a model where the optical depth profile increases from 30 to 117 au, and then decreases out to 200 au. The right set of 
panels shows a model where the optical depth increases from 30 to 70 au, as is then flat out to 200 au. 


centred on 55 and 130 au, with face-on optical depths of 
4.5 x 10 -4 and 1.3 x 10 -3 . The observed-model residuals in the 
right panels of Fig. [6| show that this model is a good fit to the 
PACS observations, with the only significant remaining struc¬ 
ture attributable to the high background level. The model also 


reproduces the disc spectrum well, though a lack of far-IR pho¬ 
tometry means that this goal is easily met. 

We also considered two similar extended models for the 
disc. These consist of two disc zones, that extend from n n to 
t mid? and then from r m id to r ou t- These two zones have inde¬ 
pendent power-law optical depth profiles r oc r a \ where i is 
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1 or 2 for the inner and outer zones. The absolute level of the 
zones is set such that the join at r m id is smooth. The disc in¬ 
clination and position angles are as in the two ring model, and 
allowing them to vary does not improve the model fit. 

For a first model, we set a 1 = 1 and a 2 = —1, so the 
optical depth increases from some inner edge to a peak value, 
and then decreases to the outer edge. For these assumptions we 
find n n = 30, r m id = 117, and r ou t = 200 au. To explore 
the parameter uncertainties we construct a second model with 
ai — 1.7 and a 2 — 0, for which we find n n — 30, r m id — 
70, and r ou t = 200 au. These two models are shown in Fig. 
[TJ which both have similarly good residual images, indicating 
that neither of these models should be preferred relative to each 
other, or relative to the two ring model described above. This 
two zone model clearly has some degree of degeneracy, with 
v mid shifting to account for changes in the optical depth profiles. 
The inner and outer radii are common however, suggesting that 
for continuous models the inner edge is around 30 au, and the 
outer edge around 200 au. 

In comparison to the two ring model, the two zone model 
extends to both smaller and larger radii. This is to be expected 
however; a model that was continuous only between the two 
radii found by the two ring model would have too much emis¬ 
sion from the region between the ring locations, so this concen¬ 
tration is countered by making the continuous model cover a 
wider range of radii. 

We briefly compare the prediction of our model for the 
scattered light brightness with limits set by iKrist et al.1 (2012) 
using the Space Telescope Imaging Spectrograph (STIS) instru¬ 
ment onboard the Hubble Space Telescope (HST). The limits 
decrease with increasing radius, being about 21 mag arcsec -2 
at 3.5" (55 au) and 23.3 mag arcsec -2 at 8.5 " (130 au). If 
we assume isotropic scattering and an albedo of 0.1, the in¬ 
ner component of our two belt model is too faint to detect (22 
mag arcsec -2 ), and the outer component is near the limit of de¬ 
tectability. Therefore the upper limit on the albedo is about 0.5 
for the inner belt, and about 0.1 for the outer belt. For the con¬ 
tinuous model the surface brightnesses are about one magnitude 
fainter for an albedo of 0.1, and the inner regions of such a disc 
cannot be detected even with an albedo of 0.9, at which point the 
outer disc is at the limits of detectability. Low albedoes are gen¬ 
erally the norm for s cattered light observations of debris discs 
(e.g. IKrist et alJEoioh . with rel atively strong forward scatter¬ 
ing be ing a possible explanation (iMin et alJ20ld : lMulders et ahl 
l2013h . Therefore, the STIS observations do not set sufficiently 
strong limits that the dust albedo in the HD 38858 disc appears 
any different to discs that were detected in scattered light. 

In summary, we find that several different models can re¬ 
produce the structure in the HD 38858 disc. It is not clear from 
the PACS data whether the disc is split into two separate compo¬ 
nents at about 55 and 130 au, or whether the disc covers a wide 
region from about 30 to 200 au. For the latter models, the opti¬ 
cal depth profile is not clear because it is possible to construct a 
two zone model with differing profiles in each zone, depending 
on the radial location of the join between these zones. 


4.1.3 SCUBA-2 interpretation 

As noted in section l3?L2l our SCUBA-2 images appear to suf¬ 
fer from some background contamination, most likely from the 
Orion Complex. As suggested by the contours in Fig. [2) the 850 



Figure 8. Residual 850 jum SCUBA-2 image (contours), with residual 
160 gm image (colours) in the background. The contours are at 2, 3 and 
4cr (solid are +ve, dashed are -ve). The star location is indicated by the 
white cross. The SCUBA-2 850 pm beam (13") is shown at the lower 
right. 


pm image likely also contains flux originating in HD 38858’s 
debris disc. Having derived possible disc structures using the 
PACS data, we can therefore assume that the structure is similar 
at 850 pm and attempt to verify this interpretation. 

Using the models of the disc structure derived above 
(which all yield similar results), Fig. [8] shows 850 pm residual 
contours after subtracting a model with a total disc flux of 7.5 
mJy. The model was created in the same way as the PACS mod¬ 
els, but at a wavelength of 850 pm and using a 13" Gaussian for 
the PSF. The residuals show a “stripe” of emission to the East of 
the stellar position, which extends between two regions where 
significant residual emission remains in the 160 pm image (see 
also Figs. 0 and [TJ. The choice of total disc flux is somewhat 
subjective, and was chosen so that the stripe is fairly uniform 
in brightness from North to South; we assign an uncertainty of 
2 mJy (i.e. 3x higher than the point source uncertainty). That 
we can successfully subtract a symmetric disc model centred on 
the expected stellar position, and that the residual flux has good 
correspondence with the 160 pm residual emission, shows that 
our interpretation is reasonable. It may of course be that the 
disc is asymmetric and some of the residual 850 pm flux to the 
East actually belongs to the disc. However, given that the PACS 
images, that have much higher S/N, are consistent with being 
symmetric, we consider this possibility unlikely. 


4.1.4 Gas mass 


To estimate the upper limit on the total CO gas mass from 
our limit on the integrated J=2-l flux (section [3.1. 3 L we apply 
the CO excitation model described in iMatra et al. (20151) . The 
model solves the statistical equilibrium of CO rotational levels, 
taking into account both collisional and radiative excitation in 
a full non-local thermodynamic equilibrium (NLTE) approach. 
The radiation field impinging on CO molecules at the frequency 
of the transition J=2-l (230.538 GHz) is dominated by the cos- 
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mic microwave background (CMB), with negligible contribu¬ 
tion from both the star (beyond about 1 au) and from any of the 
disc models described above (regardless of location). That is, 
the radiative excitation of CO molecules is largely independent 
of their radial location. Additional excitation of the molecules 
by other means, most likely collisions with electrons freed from 
carbon ionisation following CO dissociation, will only lead to 
more J=2-l emission from a fixed mass of CO. We can therefore 
set a hard upper limit on the total CO mass around HD 38858 
of 1.5 x 1(T 4 M®. 

The age of HD 38858 is well beyond when gas of any 
kind has been detected around young st ars JZuckerman et al.1 
ll995l : lMoor et al.lEoill : lMatra et alJ l2015). so little is expected 


to exist in this system, and any gas that is present likely orig¬ 
inates in icy planetesimals w ithin the observed debris disc 
(e.g. IZuckerman & Sond 12012). If we assume such a second- 
generation scenario where the gas is composed of cometary 
molecules such as H 2 O, CO and their photodissociation prod¬ 
ucts, electrons (produced via CO photodissociation followed by 
C ionisation) are the main colliders and we can set stricter lim¬ 
its on the total CO mass. Such a constraint can be illustrated by 
assuming a v alue of 100 for the C/CO ratio (as observed in /3 
Pictoris . [Roberge etal . 2000), and an ionisation fraction C + /C 
of 1.7 (as estimated bv lFernandez et al.l2006l . at 100 au around a 
G-type star). For HD38858, the interstellar UV field dominates 
photodissociation of CO even at the inferred inner edge of the 
dust disc (30 au), yielding a CO lifetime of ^120 years, while 
H 2 O will be affected by the stellar radiation field, dissociating in 
3 years at 30 au, and in 40 years at 200 au. For gas released from 
planetesimals with Solar Sys tem cometary CO/H 2 O abun dance 
ratios (between 0.4 and 30%. lMumma & CharnlevlEoIlk elec¬ 
trons dominate the collisional excitation of CO. We then find a 
CO mass upper limit of 3.2 x 10 -5 M®. While this limit de¬ 
pends on the above assumptions, it shows how the limit on gas 
can change with the scenario assumed. 

We then consider what this upper limit means in terms of 
production and destruction of secondary gas in the system. CO 
is mainly destroyed via photodissociation by the interstellar UV, 
with a timescale of ~120 years. Dividing the upper limit on CO 
mass by this timescale, we set an upper limit on the CO produc¬ 
tion rate of 1.6 x 10 18 kg yr -1 . On the other hand, we can es¬ 
timate the CO production rate from the mass loss rate of solids 
through collisions in the disc (10 15 to 10 16 kg yr -1 , section 
14.1.1b . where the production rate depends on the ice/rock and 
CO/H 2 O fractions. For Solar-System-type comets, we therefore 
obtain an observational upper limit on the CO production rate 
from the JCMT observations that is at least two orders of mag¬ 
nitude higher than the upper limit from the derived planetesi- 
mal destruction rate. Thus, the limit on CO gas is not stringent 
enough to test different scenarios for the production of gas by 
planetesimal collisions. 


4.2 HD 20794 

For HD 20794 we first consider the SED, and then show what 
the non-detection of excess at mid-IR wavelengths implies for 
the radial location of the dust detected by near-IR interferom¬ 
etry. We then turn to the images and consider the disc spatial 
structure. 



Figure 9. SED for HD 20794. Dots are fluxes and triangles 3a upper 
limits. Black symbols are measured fluxes and brown symbols are star- 
subtracted (i.e. disc) fluxes. Grey dots are disc fluxes that are consistent 
with zero, and grey triangles 3a upper limits on these fluxes. The 5480K 
stellar photosphere model is shown in blue, the 80K blackbody disc 
model in red, and the star+disc spectrum in black. 


4.2.1 SED 

We follow the same procedure described above to model the 
HD 20794 spectrum, shown in Fig. [9] The photospheric pre¬ 
diction in the far-IR is well constrained by the IRAS and 
AKARI observations. We find a stellar effective temperature 
of 5480 zb 20K, and a disc temperature of SOt^K. The disc 
temperature is very poorly constrained due to the lack of disc 
detections over a range of wavelengths and because the disc is 
fairly faint. Despite the large temperature uncertainty, the disc 
fractional luminosity is reasonably well constrained to between 
10“ 5 and 10 -6 , because disc models that fit the disc spectrum 
have similar integrated luminosities. The lower end of this range 
is only an order of magnitude above the Solar System level of 
10 -7 dVitense et al.112012!) . 

Based on the disc temperature, the blackbody disc radius 
is about 10 au, or 3" in diameter, though the uncertainty is very 
large because the blackbody radius is proportional to the in¬ 
verse square of the temperature. This size prediction suggests 
however that the disc may be resolved, as was also suggested 
by the images in Fig. |4] We now consider the location of dust 
detected with near-IR interferometry, and then turn to a more 
detailed consideration of the PACS images to test whether the 
disc is indeed resolved, and if so, whether its properties can be 
constrained further. 


4.2.2 Location of the hot dust 

A near-IR excess at 1.65 gm (H- band), with a disc/star flux 
ratio of 1.64 ± 0 .37 percent, was recently detected around 
HD 20794 bv lErtel et al.l (l2014h . Such excesses are fairly com¬ 
mon and are thought to arise from hot dus t located very near 
to the parent star (e .g. iDefrere et al.l 1201 ll : lAbsil et all 1 20131 : 
iLebreton et al.l 1201 3h . Their origin remains a mystery; dust 
should be removed by radiation forces on orbital timescales, and 
local planetesimal populations that resupply the dust would be 
rapidly depleted to undetectable levels by collisions. There is 
also no support for a non-local origin, in particular because the 
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Figure 10. Scattering and absorption properties for silicate/organic 
grains as a function of grain diameter. The dotted and dot-dash lines 
show Qabs and Qsca at H-band, and the dashed line shows Q a b s at 24 
fim. The solid line shows the albedo Q sca / (Qsca + Q a bs)- The H-band 
albedo relative to the 24 gm emission efficiency is maximised for grains 
about 0.5 fim in size. 


hot dust brightness does not correlate with other system proper¬ 
ties such as the presence o f outer dust belts (but may correlat e 
positively with stellar age. lAbsil et al.ll201~3l : lErtel et al .1120141) . 
It may be that the phenomenon is more closely related to stellar 
astrophysics, for example being very small dust trapped by the 
stellar magnetic field dSu et al.l2013t) . though in such a scenario 
the dust may still be supplied by known debris disc proces ses. 

Among the systems detected by lErtel et aD ( 2014 ). the 
spectral slopes seen across the 77-band (1.55 to 1.75 /xm) 
are relatively flat compared to the stellar flux. In the case of 
HD 20794 dust cooler than about 500 K could be ruled out, 
and while hot dust remains a possibility, the spectrum is also 
consistent with arising from scattered starlight. This latter pos¬ 
sibility means that the dust could lie anywhere within the 400 
mas PIONIER field of view, and for a given surface area of dust 
therefore extend farther from the star than if the emission is ther¬ 
mal, perhaps even into the habitable zone. No excess was seen 
at 24/xm for HD 20794, so the upper limit on the dust seen at 
this wavelength provides useful constraints on where the hot 
dust lies. Deriving this constraint however requires exploration 
of the grain properties, in particular the albedo and the emission 
or scattering efficiencies at H -band relative to the emission ef¬ 
ficiency at 24 fi m. 

The grain scattering (Q sca ) and emission/absorption 
(Qabs) efficiencies are shown in Fig. ITOl (the emission and ab¬ 
sorption efficiencies are equal). Lines show these values as a 
function of grain diameter for non-porous grains composed of 
1/3 silicates and 2/3 organics JAugereau et all 19991) at H-band 
and 24 /im. The albedo, uj = Q SC a/(Qsca + Qabs) at H-band, 
is also shown. For small grains the albedo becomes very small 
because grains absorb better than they scatter, and the maxi¬ 
mum ratio between the H-band albedo and the 24 gm emis¬ 
sion efficiency is reached for grains about 0.5 gm in size. These 
properties set limits on the grain size if the PIONIER emis¬ 
sion originates in scattered light, because no more than 100% 
of the starlight can be intercepted by the dust. That is, assuming 
isotropic scattering the fraction of starlight intercepted by the 


Temperature (K) 
587 



Figure 11. Limits on the fraction of starlight intercepted and the lo¬ 
cation of 0.5 fi m size dust grains around HD 20794 detected with PI¬ 
ONIER. The near horizontal solid line shows a locus of possible dust 
locations if the detection originates in scattered light. The dotted line 
shows a locus of dust locations if the detection originates in thermal 
emission from blackbodies. The dashed line shows the 3 a upper limit 
set by the MIPS 24 gm non-detection. The vertical stripes show the 
locations of the planets. 


dust / ca p = 0.0164 /cj < 1 so from Fig. [lO] the grains must 
be larger than about 0.1 gm for the assumed composition. As 
an example the locus of possible dust locations for grains 0.5 
gm in size is shown as the solid line in Fig. [TT] Without further 
constraints, the dust could lie anywhere along this line, the ra¬ 
dius only being restricted by the PIONIER field of view. At each 
point along this locus however, a corresponding amount of ther¬ 
mal emission is expected, so some locations can be ruled out by 
an analogous locus derived from the 24 gm non-detection. 

A constraint on the hot dust location based on the 24 gm 
non-detection is shown as the dashed line in Fi g. [TT! This line 
was computed using equation (11) from IWvat t ( 12008 3. but the 
dust temperature and Q a b s (= 1/A/0 were calculated self- 
consistently for 0.5 gm grains. The dust must lie below the 
dashed line, so if the PIONIER emission originates in light 
scattered off 0.5 gm dust with our assumed composition it lies 
within 0.03 au. For smaller dust the albedo becomes smaller 
faster than the decrease in 24 gm emission efficiency, so the 
radial constraint becomes stronger. For larger dust the 24 fim 
emission efficiency increases faster than the albedo so again the 
radial constraint is stronger. 

An additional issue for a scattered light origin is forward 
(or backward) scattering, which could be relevant here as we in¬ 
fer particle sizes similar to or larger than the PIONIER-observed 
wavelength of 1.65 gm. For a near edge-on disc the total flux 
can exceed that from isotropically scattering dust, and thus the 
solid line may be overestimated in Fig. 10. The enhancement in 
the disc flux required to alter our conclusions is large however; 
the solid line would need to move down by two orders of mag¬ 
nitude for the dust location to be allowed near 1 au. To make 
a simple estimate, we computed the integrated emission from a 
circular dust ring using the lHenvev & Greenstein ( 1941 ) phase 
function. For the emission to exceed that compared to isotropic 
scattering by an order of magnitude requires very strong for¬ 
ward scattering (g > 0.93) from a near-perfectly edge-on disc. 
In this case however, almost all of the emission originates from 
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a ten-degree arc in the ring. This kind of emission would mimic 
the off-centre point-like emission from a companion, rather than 
that of a centro-symmetric disc, and the former has already been 
ruled out for HD 20794 (tMarion et al .1120141) . This illustration is 
of course very basic, but shows that forward scattering is very 
unlikely to result in a brightness enhancement that is both con¬ 
sistent with the PIONIER observations, and results in a much 
smaller fraction of starlight intercepted by the dust compared to 
what would be inferred assuming isotropic scattering. We there¬ 
fore conclude that our constraints on the dust location (if it orig¬ 
inates in scattered light) are robust. 

The green dotted line in Fig. QT] shows the location of the 
hot dust if it originates in thermal emission, computed in the 
same way as at 24 /im. Inside about 0.5 au the line lies below 
that for scattered light, so the same signature can be created 
from thermal emission but with several times less surface area 
in dust. The 24 gm constraint still applies, so for thermal emis¬ 
sion from 0.5 gm grains the dust lies within 0.1 au. This con¬ 
straint is similar for smaller grains; Fig. [Tol shows that the ratio 
of Qabs at H-band and at 24 gm is similar for grains smaller 
than 1 gm, and thus for smaller dust both lines move upwards 
in concert. For larger grains the lines move downwards slightly 
and the radial constraint moves inwards until Q&bs = 1, and 
the constraints then do not change with increasing grain size 
(i.e. they behave as black bodies above about 10 /mi in size). 
These radial constraints on the dust location are similar for other 
compositions, for example the possibility that the grains con¬ 
tain carbon is motivated by the higher sublimation temperature 
iLebreton et al .120131) . The constraint moves inwards somewhat 
for compositions where Q a b s at 24 gm relative to H-band be¬ 
comes larger than the example here (e.g. due to a strong spectral 
feature for small silicate grains). 

If the dust is long-lived, which is suggested by the rel¬ 
atively common detection of such dust around other stars, it 
would be unlikely to reside near the planets, which would 
rapidly accrete or scatter the dust (and any parent planetesi- 
mals that could act as a local mass reservoir). The 24 gm limits 
rule out significant H-band emission beyond the planets, so the 
dust most likely lies interior to them. Thus, our analysis based 
on non-detection of excess emission with MIPS at 24 /mi sug¬ 
gests that the hot dust seen around HD 20794 with PIONIER 
lies within ~0.1 au of the star (<20 Solar radii). 


4.2.3 Resolved models 

Returning to the PACS images of HD 20794, we first test 
whether the disc is resolved by constructing a point source 
model of disc emission. The components of this model are a 
source at the star location, the NE background source, and the 
bright SE background source. As in Fig. |4j we take the motion 
of the star between the 100 and 70 /xm epochs into account. The 
point source at the star location represents the star plus an unre¬ 
solved disc. This is in fact the PSF fitting method used to derive 
the point source fluxes quoted in section 13.21 and the left three 
columns in Fig. [I2l show the results of this model. 

The left column shows the observations, the second col¬ 
umn the point source model (convolved with the PACS beam), 
and the third column the residuals after subtracting this model 
from the observations. The best-fit point source fluxes (after 
subtraction of the stellar flux) are 10.9, 12.3, and 8.9 mJy at 
70, 100 and 160 gm, where the uncertainties are the same as the 


point source uncertainties of 5, 3, and 3 mJy given in section 
13.21 In general the model is very good, finding significant fluxes 
(i.e. an excess), but notably leaves residuals on either side of the 
star at 100 gm, which is the deepest image. 

We therefore allow the disc component to be resolved, 
modelling it as a simple narrow ring of width 5 au with free 
parameters of disc radius, inclination, and position angle. The 
best fitting resolved disc model is shown in the two rightmost 
columns of Fig. [12] where the fourth column shows the disc 
model at high resolution, and the fifth column shows the resid¬ 
uals for this model. The model has a radius of 24 au, an inclina¬ 
tion of 50°, and a position angle of 8° E of N. 

This resolved model accounts for the residual flux that was 
present on either side of the star for the point source model, 
but given that this difference lies in only a few pixels that were 
about 3cr outliers, the extra complication added by a resolved 
model is only marginally warranted, and therefore uncertainties 
on the disc parameters are moot. Despite this issue, it seems 
likely that the disc is resolved with a high inclination, with a 
diameter of around 50 au (10"), and as one of the nearest Sun¬ 
like stars with a disc and planets is a high priority for future 
observations. 

Given that both 61 Vir and HD 38858 host discs with sig¬ 
nificant radial extent, we can consider what the observations 
rule out for such extended discs around HD 20794. By anal¬ 
ogy with the two-belt model for HD 38858, we added a second 
5 au wide belt at a radius of 50 au to the resolved HD 20794 
model. To remain undetected the optical depth in this belt must 
be less than about half the value in the inner belt. For HD 38858 
the outer belt has three times greater optical depth than the in¬ 
ner, and for 61 Vir the two belts have similar optical depths. A 
basic conclusion is therefore that the material in the disc around 
HD 20794 is the most centrally concentrated of the three. 


5 PLANETARY SYSTEM STRUCTURE 

We now discuss our results within the context of the wider plan¬ 
etary system structure. We begin by presenting a brief update on 
the radial velocity (RV)-discovered planet around HD 38858, 
and then RV limits on undiscovered planets around HD 20794, 
61 Vir (HD 115617), HD 69830, and HD 38858. A detailed 
analysis of the High Accuracy Radial velocity Planet Searcher 
(HARPS. iMavor etaD 12003 ) results for these systems will be 
presented in Marmier et al. (in preparation). We then go on to 
discuss possible relations between the disc structure and known 
and unknown planets around these stars. 

HD 38858b was announced by IMavor et aO (1201 lh . as a 
planet with a minimum mass of 30.55 M® orbiting at 1 au, 
using 52 RV measurements taken over eight years. With addi¬ 
tional data it has become clear that this signal was spurious and 
arises due to an alias of the magnetic cycle of the star, which 
has a period of 2930 days. The additional data (now 96 mea¬ 
surements with several years more coverage) have however re¬ 
vealed a planetary signal with a period of 198 d= 1 days (0.64 
au) and minimum mass of 12 d= 2M®. 

5.1 Limits on undiscovered planets 

To consider the structure of the planetary systems in light of 
our results, we now present the limits on undiscovered plan¬ 
ets set by HARPS data. The limits depend primarily on the RV 
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Figure 12. Disc models for HD 20794. Each row of images is a different wavelength, the top two rows are for the 70 and 160 /dm observation, and the 
bottom two rows are for the 100 and 160 gm observation. The left (first) column shows the PACS images, the second column shows the model where 
the disc is assumed to be unresolved and the third column shows the residuals for this model (i.e. first column - second). The fourth column shows a 
resolved disc model and the fifth column shows the residuals for this model. The point source model leaves residuals on either side of the star at 100 
/dm, and these are accounted for with the resolved disc model. 


precision (K in m s -1 ) and the time span over which the ob¬ 
servations were made. The limits are therefore approximately 
M p i sinz > UK M*a p \ in units of Earth masses (with M* 
in units of M© and a p i in au) for orbital periods less than the 
span of observations, and increase very steeply for longer orbital 
periods (roughly as dpi 5 ). 

Fig.[T3lshows the limits on unseen planets derived from an 
analysis of the HARPS RV data. In each case the planetary sig¬ 
nals (and stellar signals where necessary) were first subtracted 
from the time series. Then, for each point in the mass-period 
diagram, a corresponding circular orbit is added to the residu¬ 
als, with 20 phase values (spaced equidistantly between 0 and 
2tt). For each phase, this process is repeated 1000 times with 
the residuals randomly permuted, and we consider that the syn¬ 
thetic planet is detected only if we see a signal in the resulting 
periodogram that is higher than the 1 per cent false alarm prob¬ 
ability limit at all phases. 

All four systems have been observed for a time span of 
about 10 years, which is reflected in the planet sensitivity in 
terms of semi-major axes. The differences in mass sensitivity 
( K ) are almost entirely due to differences in the number of mea¬ 
surements, with 96 for HD 38858, 239 for HD 69830, 209 for 



0.1 1.0 10.0 
Semi-major axis 


Figure 13. Planet minimum masses and semi-major axes for HD 20794, 
HD 69830, 61 Vir, and HD 38858 (dots, joined by dot-dashed, dashed, 
and dotted lines respectively, with no line for HD 38858b). Solid lines 
show the 3cr limits on planets set by the HARPS RV data (the legend 
has the same vertical ordering as the lines). 
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61 Vir, and 461 for HD 20794. Fig. I~i~3l shows that planets in- 
terior to HD 38858b with masses greater than about 10 M® 
could have been detected in this system. Thus, while 61 Vir 
or HD 69830-like planetary systems can be ruled out around 
HD 38858, lower mass planets (as seen around HD 20794) may 
still be present in this system. In all cases Saturn-mass planets 
within about 10 au could have been detected. 


5.2 Radial structure and history 


Combining the RV limits, our above disc modelling results, 
those for 61 Vir from Wvatt et°aD ( 20121) . and those for 
HD 69830 from lLisse et al.1 ( 20071) . Fig. fl4l shows cartoon de¬ 
pictions of the radial structure in each system. As in Fig. |T3j 
planet masses and the limits on undetected planets are shown. 
The disc extents are derived from the Herschel images, or con¬ 
strained by mid-IR imaging an d modelling of the disc spec¬ 
trum in the case of HD 69830 (iLisse et al.ll2007l : ISmith et all 
2009). The hatched regions cover the approximate radial loca¬ 
tions where the dust is seen, with the solid regions for 61 Vir and 
HD 38858 showing the location of the two-belt models. We have 
not included the hot dust seen around HD 20794 fertel et al.1 
120141) as the origin of this phenomenon is unknown, and may 
be more related to stellar physics than planetary systems. 

Contrasting the four systems, most notable is that 
HD 69830 does not possess a significant level of cool dust 
iMarshall et al.ll2014l) . suggesting that this system evolved in 
a different way to the other three, and that the dust at 1 
au is perhaps short-lived and directly linked to one or more 
of the three planets. Th is system has been d iscussed ex- 
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tensively elsewhere (e.g. [B eichman et al.1 
20061: Llsse et al.| 12007 : Wvatt et al.1 12007 al : iPavne et al.ll2009l : 
Smith et al.ll2009( : lBeichman et al.ll201ll) . so here we focus on 


the remaining three. 

Our main goal is to explore links between the planets and 
the disc structure, whether these links provide potentially use¬ 
ful information about the system histories, and which observa¬ 
tions are needed to distinguish among different possibilities. We 
first highlight some relevant similarities and differences among 
these three systems: The discs around HD 20794, 61 Vir, and 
HD 38858 have similar inner radii of around 20-30 au, and (as¬ 
suming that the minimum masses are similar to the true planet 
masses) the planets around 61 Vir and HD 38858 are more mas¬ 
sive than those around HD 20794. The discs around 61 Vir and 
HD 38858 extend to large radii (>100 au), and while the con¬ 
straints are much poorer, the above modelling shows that the 
HD 20794 disc is at least relatively depleted outside ~30 au. 
From the SEDs and images of 61 Vir and HD 38858 we can 
also be confident that the regions inside ~30 au are significantly 
depleted of debris compared to the levels seen between 30-200 
au. Being very faint and marginally resolved at best, the relative 
constraints for HD 20794 on the regions inside the detected belt 
at ~20 au are poor. 

Ideally we would like to link the properties of the 
known planets with the observed disc structure, thus avoid¬ 
ing the invocation of unseen planets with uncertain proper¬ 
ties. Secular perturbations provide one way that planets can 
reach over large radial_distances and affect planetesimal belts 
(e.g. lMustill & Wvattll2009l) . However, the secular perturbation 
timescale for the outermost planet to affect the disc inner edge 



Radius (au) 

Figure 14. Cartoon depictions of the four nearby G-type systems with 
low-mass planets and debris disc detections. Increased planet masses 
around HD 38858 and 61 Vir assuming coplanarity with the resolved 
disc are included (open circles), but not visible for 61 Vir due to the 
small correction. RV limits from Fig.[l3]are shown as solid lines. The 
regions occupied by the discs are shown to the right. The uncertainty in 
the structure is indicated for HD 38858 and 61 Vir by including the two- 
belt model (solid region) and the extended models (hatched region). For 
HD 69830 and HD 20794 the approximate location is indicated by the 
hatched region. 
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is >10 Gyr for all three systems, so such a link is not possible 
given the current orbits of the known planets. 

In the absence of strong effects from secular perturba¬ 
tions, we can still make a speculative link between the masses 
of the known planets and the radial width of the discs, which 
are larger for the two discs in systems with more massive 
planets. Super-Earths may have difficulty forming in situ (e.g . 

Zhou et al.l2005l:lAlibert et al.l2006l : [Raymond & Cossoull20l4 

Inamdar & Schlichtim] 20l4) . so could have formed farther out 

where the discs are observed (i.e. near 20 au), imprinted some 
structure on the disc, and then migrated inwards to their ob¬ 
served locations. In this kind of scenario the planet mass de¬ 
pendence of the disc radial width arises because more massive 
planets excite planetesimals to greater eccentricities. That is, at 
least the outer planets around HD 20794, 61 Vir, and HD 38858 
all formed near the observed inner disc edges, but the >10 M 0 
outer planets around 61 Vir and HD 38858 were better able to 
scatter planetesimals onto eccentric orbits compared to the <10 
M 0 planets around HD 20794. The depletion of disc material 
between the current planet and disc locations is explained in this 
scenario by i) the use of this material to form the observed plan¬ 
ets, and/or ii) the shepherding of this mater ial inw ards as the 
super-Earths migrated towards the star (e.g. Zhou e t~aPl2005l : 

iKennedv & KenvorJl2008l: llzidoro etakll 20141) . 

By analogy with the Solar System, the radial width of the 
discs around 61 Vir and HD 38858 would therefore be inter¬ 
preted as arising from “scattered disc”-like populations of plan¬ 
etesimals on eccentric orbits. A prediction of this scenario is 
that the radial width of the disc as observed in small dust with 
Herschel is a property inherited from the larger planetesimals, 
and therefore that the disc extent should also appear to be wide 
when larger grains that are less affected by radiation forces are 
imaged at mm wavelengths. This appearance would be in con¬ 
trast to discs such as Pic and AU Mic, where the planetesimals 
occupy a relatively narrow “birth ring” and the small grains that 
are strongly affected by radiation forces reside in a wider halo 
dStrubbe & Chian g||2006h . 

While this scenario is convenient in that it does not require 
undetected planets at larger semi-major axes, it is by no means 
unique and the sensitivity limits in Figs. IT3l and fl4l show that 
the constraints allow >10-20 M 0 planets beyond about 5 au, 
and >Saturn-mass planets beyond 10 au0 Given the total mass 
of solids in a “typical” planet-forming disc is roughly 10-100 
M 0 , and the gas mass ^>100 tim es greater (e.g. lWeidenschillina 
1 19771 : 1 Andrews & Williamsll2005h . the existence of more plan¬ 
ets is certainly possible given the mass budget. Another possi¬ 
ble scenario is therefore that the known planets formed closer 
to the sta r, or at least from materia l that originated closer to the 
star (e.g. lHansen & Murravl 12012k and that extra planets may 
(or may not) exist in the regions where the RV limits are poor. 
In this case any links between the masses of the known super- 
Earths and the disc structure would be far more tenuous, and 
likely non-existent. An obvious question is then whether unde¬ 
tectable planets are responsible for the different radial widths 
and the disc depletion inside 20 au, and which signatures can 


3 Being old stars, these systems are not high-priority targets for high- 
contrast direct imaging surveys because the sensitivity to ~Gyr-old gi¬ 
ant planets is poor; only HD 69830 has been observed recently, and was 
targeted based on the presence of the disc tJanson et~ahll2013l) . 


distinguish planet formation and/or planet clearing from other 
depletion scenarios. 

Firstly, planets are not necessarily needed to clear the inner 
disc regions. Most simply it may be that planetesimal formation 
was less efficient in the 5-20 au region than elsewhere and that 
these regions are naturally empty. In the event that planetesi¬ 
mals did form in this region, IWvatt et all d2012h show that the 
disc around 61 Vir may simply be collisionally depleted from 
the inside out, and given the uncertain stellar ages, the same 
argument can be applied to HD 20794 and HD 38858. This sce¬ 
nario predic ts a characteristic slope for the brightnes s of the disc 
inner edge (IWvattl 120081 : IKennedv & Wvattll2010h . The main 
prediction if the inner disc edges are instead truncated by un¬ 
seen planets is t hat the edge should be sharper tha n for colli- 
sional depletion JOuillenl2006l : lMustill & Wvattl2012h . and that 
any relation between the known planet masses and the disc ex¬ 
tent is coincidental. The best way to tell how the inner disc re¬ 
gions were depleted is therefore to resolve the inner disc edge, 
a task best suited to high resolution optical or interferometric 
mm imaging (e.g. lKalas et al.ll2005l : lBolev et al.ll2012h . The ra¬ 
dial disc extent may of course also be related to the masses of 
these unseen planets, but unless these planets are detected or 
inferred to exist there is little power to tell whether their influ¬ 
ence on the width was significant. However, their possible pres¬ 
ence highlights other scenarios. The scattered disc component 
of the Kuiper belt formed as planet esimals were scattered b y 
the outward migration of Neptune (iPuncan & Levisonlll997h . 
providing an alternative explanation for the wider extent of the 
discs around 61 Vir and HD 38858 noted above, but no reason 
(aside from lower planet masses or the lack of such a process) 
for the narrower belt around HD 20794. The additional planets 
required in this scenario would be very difficult to detect, so a 
key signature that would distinguish among these scenarios is 
that outward mig ration captures planetesimals in mean motion 
resonances JWvattll2003h . with the resultant clumpy structures 
expected to be detectable in high resolution sub-mm images 
feitel et al.l2012h . 

In summary, the greater disc radial widths and planet 
masses around 61 Vir and HD 38858 compared to HD 20794 are 
suggestive of a link between the two; the wider extent may be 
a signature of a scattered disc population of planetesimals that 
was excited by the planets before they migrated to their current 
locations. Being based on three systems this scenario is clearly 
speculative and other scenarios are of course also possible, so 
such ideas provide general motivation for further characterisa¬ 
tion of these and other systems. In particular, high resolution 
observations that measure the shape of the inner disc edges and 
any azimuthal asymmetries will provide the best constraints on 
the structure and origins of these systems. 


6 CONCLUSIONS 

We present Herschel and JCMT observations of the debris 
discs in the Sun-like planet hosting systems HD 38858 and 
HD 20794. The HD 38858 disc is well resolved and spans a 
wide range of radii, from approximately 30-50 au to 130-200 au 
depending on whether the dust lies in two narrow rings or a wide 
continuous belt. We present a probable 850 /am imaging detec¬ 
tion of the disc and set an upper limit on the level of CO J=2-l 
emission. The HD 20794 disc appears robust photometrically, 
but is marginally resolved at best. The disc has a probable radius 
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of a few tens of au and is among the faintest known, with dust 
levels roughly one order of magnitude above the Edgeworth- 
Kuiper belt in the Solar System (e.g. lVitense et al.ll2012r) . Using 
the non-detection of excess emission at 24 gm we show that the 
hot dust seen with near-IR interferometry probably lies interior 
to the innermost planet at 0.1 au. 

With the addition of 61 Vir and HD 69830, we consider 
the sample of four nearby G-type stars with low-mass planets 
and debris disc detections. We present new limits on undiscov¬ 
ered planets in these systems using HARPS radial velocity data, 
showing that regions beyond about 5 au may contain undiscov¬ 
ered >30 M® planets. With the exception of HD 69830, the 
detected disc components lie well beyond the known planets in 
these systems and various different scenarios can explain the 
observed system structure. A scenario that does not invoke ad¬ 
ditional planets is that the extent of the observed discs around 
HD 20794, 61 Vir, and HD 38858 is related to the mass of the 
planets, which are greater around the latter two stars. Other sce¬ 
narios are possible however, and high resolution imaging that 
can determine the shape of the disc inner edges, and the extent 
of the parent planetesimals, will help distinguish among them. 
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